Is cosmology compatible with sterile neutrinos ? 
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By combining data from cosmic microwave background (CMB) experiments (including the recent 
WMAP third year results), large scale structure (LSS) and Lyman-a forest observations, we con- 
strain the hypothesis of a fourth, sterile, massive neutrino. For the 3 massless + 1 massive neutrino 
case we bound the mass of the sterile neutrino to m s < 0.26eV (0.44eV) at 95% (99.9%) c.L. These 
results exclude at high significance the sterile neutrino hypothesis as an explanation of the LSND 
anomaly. We then generalize the analysis to account for active neutrino masses (which tightens the 
limit to m a < 0.23eV (0.42eV) and the possibility that the sterile abundance is not thermal. In 
the latter case, the contraints in the (mass, density) plane are non-trivial. For a mass of > leV or 
< 0.05eV the cosmological energy density in sterile neutrinos is always constrained to be uo v < 0.003 
at 95% c.L. However, for a sterile neutrino mass of ~ 0.25eV, cu v can be as large as 0.01. 

PACS numbers: 



Introduction. — Recent cosmological data coming 
from measurements of the Cosmic Microwave Back- 
ground (CMB) anisotropies (see e.g. 0,Q|)> galaxy 
clustering (see e.g. jjj) and Lyman-alpha forest clouds 
are in a spectacular agreement with the expec- 
tations of the so-called standard model of structure 
formation, based on primordial adiabatic inflationary 
perturbations and a cosmological constant. 

Since the model works so well, the ambitious idea 
of using cosmology to test aspects of particle physics 
is becoming a reality. An excellent example of this 
comes from the new cosmological constraints on neu- 
trino physics. 

Cosmological neutrinos have a profound impact on 
cosmology since they change the expansion history of 
the universe and affect the growth of perturbations 
(see la for a detailed account). Recent analyses (see 
e.g. lit*]) have indeed constrained the neutrino mass 
in the context of three-flavour mixing to be m„ < 
0.16eV (m u < 0.45eV without Lyman-a forest data) 
with a greater accuracy than laboratory beta decay 
experiments which suggest m v < 2.2eV (see @ and 
references therein). 

A possible discrepancy between cosmology and beta 
decay or neutrino oscillation experiments might pro- 
vide valuable information for the presence of system- 
atics or new physics. At the moment, the claimed 
and highly debated detection of a neutrino mass in 
the range 0.17eV < < 2.0eV at 99% c.L from 
the Heidelberg-Moscow double beta decay experiment 
is at odds with the cosmological bound. 

While the neutrino masses are very difficult to mea- 
sure experimentally, mass differences between neu- 
trino mass eigenstates (mi ,7712,7713) have now been 



measured in oscillation experiments. Observations of 
atmospheric neutrinos suggest a squared mass differ- 
ence of Am 2 ~ 3 x 10 _3 eV 2 , while solar neutrino 
observations, together with results from the Kam- 
LAND reactor neutrino experiment, point towards 
Am 2 ~ 5 x 10 _5 eV 2 . The two measured mass dif- 
ferences are easily accomodated in simple extensions 
of the Standard Model by giving masses to at least 
two of the neutrinos. If these masses are greater than 
~ 0.1 eV, all three neutrinos must be nearly degener- 
ate, with small differences accounting for the observa- 
tions. 

Results from the Liquid Scintillator Neutrino De- 
tector (LSND) || challenge the simplicity of this pic- 
ture. The LSND experiment reported a signal for 
Pfj, — ► v e oscillations in the appearance of v e in an 
originally beam. To reconcile the LSND anomaly 
with results on neutrino mixing and masses from at- 
mospheric and solar neutrino oscillation experiments, 
one needs additional mass eigenstates. One possi- 
bility is that these additional states are related to 
right-handed neutrinos, for which bare mass terms 
(MvrVr) are allowed by all symmetries. These would 
are sterile, i.e. not present in SU(2)l x U(l)~ in- 
teractions. The "3 + 1 sterile" neutrino explanation 
assumes that the — > v e oscillation goes through 
Vfi — > v s — > v e . The additional sterile state is sep- 
arated by the three active states by a mass scale in 
the range of 0.6eV 2 < Am 2 SND < 2eV 2 . Constraints 
from long baseline experiments are threatening this 
interpretation H E lH [H El Q ; it is possible 
that more than one sterile neutrino is necessary to ex- 
plain LSND The LSND signal will be soon tested 
by the MINI-BOONE experiment, expected to release 
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the first results at the beginning of the next year. 

In the meantime, given the increased quality in the 
data, it is timely to test the sterile neutrino hypoth- 
esis using cosmological observations. Several recent 
analyses have already provided interesting cosmologi- 
cal constraints on a fourth massive neutrino |4i ll6lll7l |. 
Here we generalize these in several ways: first, while 
previous work has considered the case of 3 (massless) 
+ 1 (massive) neutrino, here we also allow masses for 
the standard 3 neutrinos, as required by oscillation 
experiments. Second, we use updated cosmological 
datasets,including the new results from the WMAP 
satellite 0,13 and BOOMERAG-2K2 experiment GH. 
Finally, the strength of the interactions of a neutrino 
determines its cosmological abundance. Given how 
little we know about sterile neutrino interactions (or 
mass mixing), it therefore seems reasonable to allow 
the sterile abundance to be a free parameter. Of 
course, if a sterile neutrino can have any abundance 
(including zero), there is no mass limit. However, we 
will see that the constraints in the (mass, density) 
plane are highly non-trivial. 

Cosmology. — The three active neutrinos interact 
via the well- measured weak interactions. These inter- 
actions ensure that they were in thermal equilibrium 
at early times until they decouple from the primordial 
plasma slightly before electron-positron annihilation. 
After decouple, they maintain an equilibrium distri- 
bution of a massless fermion with a temperature lower 
than the photon temperature by a factor of (4/11) 1 / 3 . 
This introduces a well-known relation between the en- 
ergy density of the active neutrinos and their total 
mass: 
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where p cr is the critical energy density, h parametrizes 
the Hubble constant via Hq = lQOh km sec -1 Mpc , 
and here and throughout m v refers to the sum of all 
active neutrino masses. So, for example, if the three 
neutrinos are nearly dengerate, they each have a mass 
approximately equal to m„/3. 

While sterile neutrinos, by defintion, do not have 
weak interactions, they are not pure mass eigenstates. 
As such, oscillations in the early universe can trans- 
form the thermal active neutrinos into a sterile neu- 
trino 0, H3] . Thermalization occurs if 



ZWlsnd sin 4 9 > 3 x 10~ 6 eV 2 
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where 9 is an effective mixing angle. In the simplest 
models with one sterile neutrino, this condition is sat- 
isfied, so u s — 0.0106(m s /eV), but there are many 
ways of evading thermalization |2lj |. Indeed, if one 
light sterile neutrino exists, there is every reason to 
expect one or two more and these considerably com- 
plicate the thermalization analysis. It is, for example, 
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FIG. 1: The epoch of equality a oq as a function of mass of 
sterile neutrino and its energy density. The non-relativistic 
matter density here is fixed to ui m = 0.15, so that in the 
standard 3-neutrino model, oeq = 2.82 x 10 -4 . Notice 
that, at fixed u> 3 , a eq rises very rapidly for lower masses 
since the neutrinos behave as radiation. Thermalized neu- 
trinos lie along the dashed curve. 



possible to have super-thermal abundances if a heavier 
sterile state decays at relatively late times to a lighter 
state. In short, if a sterile neutrino exists, its cosmo- 
logical density is much more uncertain than that of 
the active neutrinos. 

Sterile neutrinos influence the development of inho- 
mogeneities and anisotropies in the universe by chang- 
ing the epoch of equality and by suppressing pertur- 
bations via freestreaming. The epoch at which the en- 
ergy density in non-relativistic matter equals that in 
radiation dictates when structure begins to grow. This 
leaves an imprint on the matter power spectrum [22^ : 
there is a peak at the scale equal to the horizon at 
the epoch of equality. If this epoch is close to recom- 
bination, the residual radiation causes gravitational 
potentials to decay, and this time variation produces 
an early Intergated Sachs- Wolfe (ISW) effect, boost- 
ing the power on scales near the horizon. The main 
effect of freestreaming is a suppression of power on 
scales smaller (wavenumber k larger) than 
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with suppression proportional to u> s /Lu m , where u> m = 
Q m h 2 , and Q m is the total energy density of non- 
relativistic matter (baryons plus cold dark matter) in 
units of the critical density. 

In the standard cosmology, with three massless neu- 
trinos, the scale factor at equality is oeq = 2.82 x 
10 _4 (0.15/w m ) . A sterile neutrino is relativistic until 
its temperature drops beneath its mass, so masses of 
order an eV raise the question: what does it count as, 
matter or radiation? Since the Hubble rate scales as 
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a~ 2 in a radiation dominated universe and a~ 3 / 2 in 
a matter dominated universe, we define the epoch of 
equality as the moment when 
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This definition agrees well with the standard defini- 
tion for massless neutrinos. The dependence of a eq on 
the sterile neutrino parameters m s and uj s is plotted 
in Figure ^ This figure suggests that in the limit of 
very small m s , any appreciable lu s will be excluded 
because neutrinos behave essentially as radiation and 
shift the redshift of matter-radiation equality signifi- 
cantly, producing an unaccetably large ISW effect. 

The amount of suppression due to freestreaming in- 
creases as the density increases (from top to bottom 
in Fig.nj, but the large scales (from which constraints 
derive) cease to be affected as the neutrino mass in- 
creases (from left to right). Therefore, at fixed oj s , 
constraints from freestreaming are tighter for small 
neutrino masses. Note that this differs from the ther- 
mal case (dashed curve in Fig. GJ. In that case, the 
neutrino density increases with the mass, so there is 
more suppression at high masses. 

Data analysis and Results- To obtain contraints 
on sterile neutrino parameters, we use the publicly 
available Markov Chain Monte Carlo (MCMC) pack- 
age CosmoMC 23] . The linear perturbations engine 
CAMB [24| of the software has been generalized in 
several ways. First, we allow for a non-thermal sterile 
neutrino density. Second, we allow for the possibility 
that the active neutrinos have mass different than the 
sterile neutrino. 

In the MCMC, we sample the following 8 dimen- 
sional set of cosmological parameters, adopting flat 
priors on them: The log mass of thermal sterile neutri- 
nos logm s and w„, the energy density of 3 degenerate 
standard massive neutrinos u> v — m„/(94.1eV), the 
physical baryon and CDM densities, ojb — Vt^h 2 and 
u) c = £l c h 2 , the ratio of the sound horizon to the an- 
gular diameter distance at decoupling, S , the scalar 
spectral index and the overall normalisation of the 
spectrum, n s and A 8 , and, finally, the optical depth 
to reionisation, r r . We consider purely adiabatic ini- 
tial conditions, impose flatness, and do not include 
gravitational waves. 

We include the WMAP three-year data 0,12 (tem- 
perature and polarisation) with the routine for corn- 
rating the likelihood supplied by the WMAP team 
i3, as well as the CBI VSA 0, ACBAR 

2M and BOOMERANG-2k2 measurements of the 
CMB on scales smaller than those sampled by WMAP. 
In addition to the CMB data, we also consider the con- 
straints on the real-space power spectrum of galaxies 
from the SLOAN galaxy redshift survey (SDSS) [2^ 
and the 2dF galaxy redshift survey |3Jj and Lyman- 
alpha forest clouds [3ll l32| from the SDSS, the gold 
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FIG. 2: Effect of extra sterile neutrino on the CMB (top) 
and LSS (bottom) power spectra. Thin lines correspond 
to standard model, sterile neutrino of mass m = leV 
(dashed) m = 0.3eV (dot-dashed) and fixed sterile den- 
sity Us — 0.01. These curves are normalised to large scale 
Ci. Thick dashed and dot-dashed curves correspond to 
models, which in addition to having sterile mass have had 
dark matter density increased to match standard a cq and 
h increased to match CMB peak positions and were nor- 
malised at the first peak. Dotted vertical lines on the 
bottom plot enclose the area where LSS experiments are 
currently sensitive to with thick line normalisations cho- 
sen to illustrate the fact that the leV model is a poorer 
fit than 0.3eV model. See text for discussion. 



sample of the recent superno va ty pe la data [33] , the 
latest SNLS supernovae data |3J] and the constraints 
from the baryonic acoustic oscillations detected in the 
Luminous Red Galaxies sample of the SDSS [13 1 

The details of the analysis are the same as those in 
|3rij ] and the reader is invited to check that paper to 
examine what constraints the above datasets give for 
other models including standard 3 degenerate massive 
neutrinos case. 

If the active neutrino masses arc fixed to zero and 



1 There is a negligible overlap between the constraints from the 
2dFGRS, SDSS and SDSS LRG analysis, as there are galaxies 
in common in all three data sets. 
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FIG. 3: 1,2-a constraints on the sterile neutrino mass and 
abundance. 



the sterile neutrino abundance is thermal (similar to 
the assumptions imposed in Ref. Q), the upper limit 
on the sterile neutrino mass is 0.26eV (0.44eV) (all 
at 95% (99.9%) c.l.). Of course the active neutrino 
masses are not zero. Taking them as a free parameter 
leads to an upper limit on the sterile neutrino mass of 
0.23eV (0.42eV). This is marginally tighter than the 
m v = constraint because the limit is really on the 
sum of all neutrino masses. Fixing the active masses 
to zero allows the maximum m s . Relaxing this restric- 
tion leaves less room for a large m s . We have found 
some sensitivity to the mass difference of the sterile 
and active states (and this might be measurable with 
future data), but current data really constrain only 
the sum of all neutrino masses. 

We now generalize further and allow the sterile neu- 
trino abundance uj s to vary. Fig. shows the con- 
straints in the uj s -m s plane. Note the distinct peak 
around the region of m s ~ 0.25eV, presenting an 
allowed region of parameter space with anomalously 
large values of uj s . To the left of this peak, oeq is very 
large and the resulting ISW effect precludes agreement 
with CMB data. When m s is in the allowed regime, 
oeq would still be too large were uj m fixed. How- 
ever, a model with larger w m (~ 0.18) leads to an even 
smaller, acceptable oeq- Fortuitously, the enhanced 
cold matter density also mitigates the freestreaming 
suppression (which scales as u>c )■ At larger neu- 
trino mass (~ leV), additional cold matter would 
make <2eq too small, so u> m must be closer to 0.13 
and the free-streaming suppression becomes relevant 
again, preventing agreement with large scale struc- 
ture. This is illustrated in Fig. Here we show 
the angular CMB anisotropy and matter power spec- 
trum for different masses at fixed lu s . The suppres- 
sion due to free-streaming is evident in the power 



spectrum and clearly becomes more severe for smaller 
masses. However, increasing dark matter density to 
match the epoch of matter-radiation equality opposes 
this effect. Crucial to this interpretation is the real- 
ization that the matter-radiation equality is very thor- 
oughly measured by the present-day experiments with 
little model-dependence. The constraint can be sum- 
marised in a cq ~ (2.95 ± 0.13) x 10~ 4 . 

Conclusions- By combining data from cosmic mi- 
crowave background experiments, galaxy clustering 
and Ly-alpha forest observations we have constrained 
the hypothesis of a fourth, sterile, massive neutrino, as 
an explanation of the LSND anomaly. For the 3 mass- 
less + 1 massive thermal neutrino case we bound the 
mass of the sterile neutrino to m v < 0.26eV(0.44eV) 
at 95 (99.9) % c.l. Marginalizing over active neutrino 
masses improves the limit to m v < 0.23eV(0.42eV). 
These limits are incompatible at more than 3er with 
the LSND result 0.6eV 2 < Am^ SND < 2eV 2 (95% 
C.L.). Moreover, our analysis renders the LSND 
anomaly incompatible at high significance with a de- 
generate active neutrino scenario and viceversa. If we 
allow for the possibility of a non-thermal sterile neu- 
trino, we find that the upper limit of allowed energy 
density in the sterile neutrino is a strong function of 
mass. In particular, for m s < leV or > 0.05eV the 
cosmological energy density in sterile neutrinos is al- 
ways constrained to be oj s < 0.003, but that for ster- 
ile neutrino mass of ~ 0.25eV, ui s can be as large as 
O.OleV. 

The results presented in this paper rely on the as- 
sumption that systematics in the public datasets we 
analyzed (WMAP, Lyman-a, etc.) are under control. 
We argue that this is likely: the datasets are large 
enough that detailed systematics checks - e.g. di- 
viding the data into multiple subsets, constructing 
quiet channels that should see nothing, and cross- 
correlating different bands to reduce noise - have been 
performed. We also checked that if we drop either 
small scale CMB, LSS or Lyman-a dataset from the 
analysis, the constraints simply weaken without any 
systematic change in the results. 

The results presented here also rely on the assump- 
tion of a theoretical cosmological model based on a 
large but limited set of parameters. Extensions of 
the parameter space - e.g., inclusion of isocurvature 
modes, gravity waves or a different parametrization 
for the dark energy component - may modify our con- 
clusions. Those modifications are however not needed 
by current data and some of them may well lead to 
stronger limits. Indeed the simple cosmological model 
with only a handful of parameters does an excellent 
job explaining a wide variety of data. If the LSND 
anomaly is confirmed by MINIBOONE, we will have 
been proved wrong, and cosmologists will need to re- 
examine the entire framework on which these very 
tight constraints rest. 
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